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Abstract 

The most general initial conditions of CMB anisotropics, compatible with the presence of 
pre-equality magnetic fields, are derived. When the plasma is composed by photons, baryons, 
electrons, CDM particles and neutrinos, the initial data of the truncated Einstein-Boltzmann hi- 
erarchy contemplate one magnetized adiabatic mode and four (magnetized) non-adiabatic modes. 
After obtaining the analytical form of the various solutions, the Einstein-Boltzmann hierarchy 
is numerically integrated for the corresponding sets of initial data. The TT, TE and EE angu- 
lar power spectra are illustrated and discussed for the magnetized generalization of the CDM- 
radiation mode, of the baryon-radiation mode and of the non-adiabatic mode of the neutrino 
sector. Mixtures of initial conditions are examined by requiring that the magnetized adiabatic 
mode dominates over the remaining non-adiabatic contributions. In the latter case, possible de- 
generacies between complementary sets of initial data might be avoided through the combined 
analysis of the TT, TE and EE angular power spectra at high multipoles (i.e. t > 1000). 



1 Pre-equality CMB initial conditions 



What are the initial conditions of CMB anisotropies? In the pivotal ACDM paradigm, the initial 
conditions are taken to be adiabatic p], EJ [3j [5] . In this case the fluctuations of the total pressure 
are proportional, prior to matter-radiation equality, to the fluctuations of the (total) energy density. 
Already in the absence of large-scale magnetic fields, the initial data of CMB anisotropies are not 
exhausted by adiabatic solution. In the non-adiabatic case the fluctuations of the total pressure arise 
because of the compositeness of the pre-equality plasma. 

Standard Boltzmann solvers include the possibility of having one (or more) non-adiabatic initial 
conditions. The generalized initial conditions of the Einstein-Boltzmann hierarch^] are hereby derived 
when large-scale magnetic fields are consistently included in the pre-equality plasma. The obtained 
results complement and extend former analyses only centered around the magnetized adiabatic mode. 

A stochastic background of large-scale magnetic field can be naturally incorporated in the physics 
of the adiabatic initial conditions. The TT angular power spectral have been analytically estimated, 
within the tight-coupling approximation, in [6], [71 [8] . The results of [TJ IB] can be used for simplified 
estimates of the TT angular power spectra, and, partially j6], of the TE and EE correlations. A full 
numerical approach is however required to confront the present j4[ [5] and forthcoming [9] experimental 
data. This problem was successfully addressed and solved in [TQl EH] • Semi-analytical results (obtained 
via a different treatment of recombination and diffusive effects) agree with the full numerical calcu- 
lation: the shape of the correlated distortions of the first three acoustic peaks is correctly captured 
by the analytical discussion even if the numerical approach is intrinsically more accurate especially 
at high multipoles. This occurrence strengthen the consistency of the numerical approach and allows 
interesting analytical cross-checks. In [TQl [TT] the only initial conditions examined were the ones as- 
sociated with the magnetized adiabatic mode. This choice is prompted by the best fit of the WMAP 
data alone as well as, for instance, by the best fits of the WMAP data combined with the large-scale 
structure data [12], [13], with the supernova data [HI [15] and with all cosmological data sets. The strict 
adiabaticity of the initial conditions will be now relaxed and non-adiabatic modes will be scrutinized 
in the presence of stochastic magnetic fields. This program is technically mandatory and physically 
relevant. 

If the initial data of the Einstein-Boltzmann hierarchy are solely non-adiabatic, the measured TE 
correlations cannot be reproduced [3] . This result is already apparent from the first 200 multipoles of 
the TE power spectra where, generically, non-adiabatic modes lead to a positive correlation while a 
predominant adiabatic mode would imply, instead, a negative cross-correlation. The position of the 
ant icor relation peak can be related, in the adiabatic case, to the position of the first Doppler peak 
of the TT power spectra (i.e. £oop — 220). The position of the first anticorrelation peak of the TE 
angular power spectrum can be estimated as £ ant i ~ 3£do P /4 — 150 to first-order in the well known 
tight-coupling expansion [16] (see also [UJ HH [19]). The experimental evidence does not exclude 
that a predominant adiabatic mode could be present in combination with sub-dominant non-adiabatic 
contributions so that the overall fit to the data may even improve [201 I2T] . 

■"■By Einstein-Boltzmann hierarchy we mean the set of kinetic equations written in curved space and supplemented 
by the contribution of the gravitational inhomogeneities obeying the perturbed Einstein equations. 

2 Following a consolidated terminology, the angular power spectra of the temperature autocorrelations will be denoted, 
with stenographic notation, by TT. In analog terms EE and TE denote, respectively, the angular power spectra of the 
polarization autocorrelations and of the temperature-polarization cross-correlations. 
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If the magnetic fields are present prior to matter-radiation equality (and, a fortiori, prior to photon 
decoupling) the physical situation changes when compared with the conventional non-adiabatic initial 
conditions. Plausible questions then arise: 



• can magnetic field be physically compatible with the various non-adiabatic initial conditions? 

• in the latter cases, what are the initial conditions to be imposed to the Einstein Boltzmann 
hierarchy? 

• how do the magnetized non-adiabatic mode affect the TT, EE and TE angular power spectra? 

• is it possible to compensate the distortions induced by the magnetized adiabatic mode with an 
appropriate non-adiabatic contribution where the magnetic field is also consistently included? 

In the aforementioned list of items, the first pair of questions implies a thorough analytical discussion 
of the initial conditions of the Einstein-Boltzmann hierarchy. The remaining two points necessarily 
demand a full numerical integration which will be performed with the code devised in [101 [U] and 
hereby extended to handle initial data which are, simultaneously, magnetized and non-adiabatic. The 
code employed for the numerical integration is based, originally, on CMBFAST [f| [22l [23] (which is, in 
turn, deeply rooted on COSMICS [231125]). 

To commence, non-adiabatic initial conditions of the Einstein-Boltzmann hierarchy will be quan- 
titatively introduced. Prior to equality, the plasma consists of a baryon-lepton component, tightly 
coupled via Coulomb scattering, and supplemented by photons, neutrinos and cold dark matter parti- 
cles (CDM in what follows). The fluctuations of the total pressure p t (denoted as S s p t ) can be written, 
in general terms, as the sum of two physically different contributions, namely, 

pi 

$sPt = C 2 st 5 s p t + 5sPnad, 4, = -j, (1.1) 

Pt 

where 5 s p t is the fluctuation of the energy density and Cg t is the total sound speed across the radiation- 
matter transition^. At the right hand side of Eq. (Il.ip . the first term (i.e. Cg t 5 s p t ) parametrizes the 
adiabatic contribution. The second term at the right hand side of Ec^ (11. ip . (i.e. S s p na ^), accounts for 
the non-adiabatic pressure fluctuations which can be also written a; 
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^sPnad = E ^<kij = 7^-7 E APj( C si - 4)*% 5 ij = "~> ( L2 ) 

ij vnpt jj ftj 

where the indices i and j are not tensor indices but denote two generic species of the pre-equality 
plasma. Furthermore, in Eq. (jl.2p . c| and are the sound speeds of two (generic) species of the 
plasma; is the fluctuation of the specific entropy computed for a given pair of species and «Sy, 
as indicated, is the relative fluctuation of ftj. By definition, <Sij = — a factor 2 (included in the 



3 In CMBFAST only the CDM-radiation and baryon-radiation modes are included. In our code all entropic modes are 
be implemented in combination with a stochastic magnetic field affecting both the initial conditions and the evolution 
equations. 

4 The background geometry is characterized by a conformally flat line element ds 2 = a 2 (T)[dT 2 — dx 2 }. The scale 
factor a(r) will be often referred to its value at matter-radiation equality, i.e. a = a/a cq . Throughout the paper the 
prime denotes a derivative with respect to the conformal time coordinate r while the overdot indicates a derivative with 
respect to the cosmic time coordinate t. The conformal time coordinate r is related to the cosmic time as a(r)dr = dt. 

°In the present notations 7i = a! /a = a' /a = aH where, by definition, H = a/a is the Hubble rate in cosmic time. 
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denominator at the right hand side of Eq. (11.21) ) avoids double counting of entropy modes when the 
sum extends over all the species present in the plasma. 

Neglecting for a moment neutrinos as well as baryons and electrons, we can imagine a plasma 
dominated by photons but with a subleading contribution of CDM particles (characterized by a con- 
centration n). In this toy example the specific entropy q ~ T'^/n will then imply S = 3<5 7 /4 — 5 C (where 
<5 7 and 5 C are, respectively, the radiation and the CDM density contrasts). For the latter estimate, it 
is useful to stress that p 7 oc while p c cx n. This sort of heuristic arguments can be phrased in fully 
gauge-invariant terms with the result that@ 

* = fa = -3 tt -a, c,^ + 3^T), d.3) 

where w\ is the barotropic index of a generic species. The entropy fluctuations of Eq. (II. 3p are 
expressed in terms of the £ and £j which are, themselves, gauge-invariant and which become, in the 
uniform curvature gauge, the density contrasts of the single species^ [26] (see also [271 EE])- The 
gauge-invariant definition of entropy fluctuations goes back to the seminal contributions of [26] and 
it is commonly employed in the generalized discussion of the evolution equations of curvature and 
entropy fluctuations (see also [29|,|30] and references therein). The quantity Q of Eq. (11.31) is expressed 
in the synchronous gauge which will be consistently used in the analytical calculations and in the 
numerical implementation of the code. In Appendix [A] the truncated Einstein-Boltzmann hierarchy is 
summarized in the language of the synchronous gauge. At early times (and in the tight-coupling limit) 
the truncated Einstein-Boltzmann hierarchy will be used to deduce the magnetized initial conditions 
for the various modes. In the pre-equality plasma we have four different species and, consequently, we 
will have that: 



c 7 = e + ^, c* = z + j, Cc = e + |, c b = e + |, (1.4) 



where 5 7 , S u , <5b and 5 C are the density contrasts of the corresponding component evaluated in the syn- 
chronous coordinate system. The electrons are tightly coupled to baryons through Coulomb scattering 
and can be treated, for the purpose of the initial conditions, as a single species with approximate 
common temperature [I I] . The initial conditions of the Boltzmann hierarchy are set after neutrino 
decoupling (i.e. well before matter-radiation equality). At this epoch, Eqs. (11.21) . (II. 3p and (11.41) imply 
that iS cr is: 

3 3 
S C7 = -3(Cc - C 7 ) = T&y ~ <*c, S cv = -3(( c - („) = -6 U - 6 C . (1.5) 

Equation (11.51) defines the CDM-radiation mode. The baryon-radiation mode causes, instead, a mis- 
match between the density contrast of the baryon-lepton fluid and the density contrasts of the (two) 
relativistic species at the corresponding epoch: 

3 3 
<S b7 = -3(Cb - Ct) = 7^7 _ $b, S hu = -3(Cb - ( v ) = j$v - S h . (1.6) 

In the ACDM choice of cosmological parameters^ lu-^q <C lumo the baryon-radiation contribution is 
suppressed as uj^q/uJcq. After neutrino decoupling the neutrino fraction R u = p u /pR is smaller than 

6 The variables Ci and Cj are directly expressed in the synchronous coordinate system (see, for instance, Eq. (jA.ljl of 
Appendix [2J ■ 

7 In equally correct terms we could also argue that Ci measures the curvature perturbations on the hypersurfaces 
where the energy density of a given species is uniform. Indeed, in the synchronous gauge £ is proportional to £ which is 
related to 1Z (see the following section) which is the curvature perturbation on comoving orthogonal hypersurfaces [26] . 

8 We shall denote, in accordance with the established conventions, ujxo — ^o^xo where hg is the indetermination on 
the Hubble rate and Qxo is the critical fraction of a given species. 
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the photon fraction _R 7 = p 7 /p R = 1 — R v : 



r 7 /4\ 4 / 3 /N,A 

R ^—r- r =8 N »(n) S0 ' 681 (f)' < L7 > 



S^ = -3(C„-C) = 7(* 7 -W- (U 



The third entropic fluctuation resides then in the neutrino-photon sector, i.e. 

3 
4 1 

A further solution of the truncated Einstein-Boltzmann hierarchy consists in initializing the lowest 
multipoles by assuming that the neutrino and baryon-photon dipoles are non- vanishing but still satisfy 
the momentum constraint of Eq. (1A.3j) . This solution is generically known as the neutrino- velocity 
mode. Finally, in the terminology of Eqs. (II. 6p and (11. 8p the adiabatic mode arises when all 
vanishes identically, which implies, in the pre-equality plasma, that C7 = C> = Ct> = Cc- 

As pointed out in the past [26] (see also [31]) the most general initial data of the Boltzmann 
hierarchy can be summarized, in the absence of a magnetized contribution, by a 5 x 5 matrix [321 [33] 
where the dominant amplitude might be the one associated with the adiabatic mode. This approach 
stimulated, in recent years, various attempts of including non-adiabatic contributions on the CMB 
initial data [20J 121] (see also [31 [351 [36] ) . 

The Boltzmann hierarchy can be initialized by a mixture of modes: one of the modes can be 
adiabatic and others non-adiabatic. Predominantly adiabatic initial conditions are the ones where the 
adiabatic amplitude is larger than the various non-adiabatic amplitudes. The strategy defined in the 
two previous sentences summarizes the bottom-up approach to initial conditions. Combinations of one 
(dominant) adiabatic mode and of other (subleading) modes are not excluded but sometimes even help 
in improving the overall fit to the cosmological observables [36J. If, on top of the adiabatic solution, 
there are two or three non-adiabatic modes then it is hard, with the present data, to infer stringent 
constraints on the non-adiabatic components. 

Why is it important to include consistently magnetic fields when setting initial conditions for the 
adiabatic and non-adiabatic modes? The answer on the nature of large-scale magnetization in the 
present Universe is still under active discussion [371 ESI EH] • The scrutiny of CMB initial conditions is 
a powerful window on the possible existence of large-scale magnetic fields in the pre-equality stage. If 
the origin of the large-scale magnetic fields is primordial (as opposed to astrophysical) it is plausible to 
expect the presence of magnetic fields in the primeval plasma also before the decoupling of radiation 
from matter. CMB anisotropies are germane to several aspect of large-scale magnetization (see, for 
instance, [101 HI] for two topical reviews on the subject). In the recent past valuable discussions of 
the role of pre-equality magnetic fields have been conducted in different frameworks. For instance 
fully covariant approaches have been used to characterize the interplay between large-scale magnetic 
fields and curvature perturbations in their relativistic regime [4*2l |4"3] (see also [39J). In other studies 
the vector and the tensor modes induced by large-scale magnetic fields have been more specifically 
addressed [SI US]. The approach pursued in this paper, directly linked to the theoretical framework of 
[3 El [8], completes and extends earlier works insofar as previous studies did not address specifically the 
calculation of the standard CMB observables for different initial conditions of the Einstein-Boltzmann 
hierarchy. The present approach wants to fill this gap and bring the analysis of magnetized CMB 
anisotropies to the same standard of the conventional case where large-scale magnetic fields are not 
included. Of course, in the past, various interesting effects have been pointed out for specific magnetic 
field configurations. Elegant formalisms have been also explored. While this was not a useless exer- 
cise, the quality of the present observational data clearly demands new and more sound theoretical 
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calculations both at the analytical and numerical level. At the same time the improved theoretical 
and numerical tools should also delicately improve on the case where magnetic fields are absent: in the 
opposite case the comparison with the putative ACDM model (and its possible improvement) would 
be much more cumbersome. 

The plan of the present paper is then be the following. In Section [2] the truncated Einstein- 
Boltzmann hierarchy is discussed in the presence of large-scale magnetic fields. In Section [3] the 
magnetized adiabatic mode is swiftly reviewed. Sections HI [5] and [6] will be devoted, respectively, to 
the magnetized CDM-radiation mode, to the magnetized baryon-radiation mode and to the entropic 
mode of the neutrino sector. In Section [7] we will discuss the case of mixtures of initial data when the 
magnetized adiabatic mode dominates over the other (non-adiabatic) contributions. Concluding dis- 
cussions are collected in Section El Without indulging in idle details, the relevant analytical tools have 
been summarized, within the synchronous gauge and within our set of conventions, in the Appendix 
lAl This choice makes the present script reasonably self-contained. 

2 Initial data of the Einstein-Boltzmann hierarchy 

The Einstein-Boltzmann hierarchy, appropriately truncated to the lowest multipoles, is initialized by 
complying with various intermediate steps: 

• the baryon-photon evolution equations have to be solved the tight-coupling approximation^]; 

• the neutrino hierarchy should be truncated to the quadrupole (or, depending upon the specific 
mode, to the octupole); 

• after solving the evolution equations for the CDM evolution, the metric fluctuations can be 
computed (for instance in the synchronous gauge). 

The momentum and Hamiltonian constraints (stemming, respectively, from the (Oi) and from the (00) 
components of the perturbed Einstein equations) have to be consistently enforced on the set of initial 
data. The Einstein-Boltzmann hierarchy must be initialized deep in the radiation-dominated stage of 
expansion and for typical wavelengths larger than the Hubble radius at the corresponding epoch. The 
two mentioned physical limits define, as we shall see in a moment, the relevant expansion parameters 
of the problem. The solution of the background Friedmann-Lemaitre equations (i.e. Eq. (1A.7j) ) across 
the radiation-matter transition (and for a spatially flat Universe) reads 

« = - = x 2 + 2x, x = l ri = (^)"m P c, (2.1) 

a eq T X H V ll M0 V0.134/ 

where a eq is the scale factor at the equality, i.e. the moment when non-relativistic matter and radiation 
contribute equally to the total energy density of the plasma. For a < 1 (i.e. a <C a eq ) the plasma is 

9 The tight-coupling expansion implies, to zeroth-order, that the photon-baryon velocity are equal (see [16] for the 
pioneering work on the semi-analytical description of scalar CMB anisotropics in the tight coupling limit). This approx- 
imation is used, at early times, also in the Boltzmann solvers to avoid numerical instabilities related to the stiffness of 
the problem. 
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dominated by radiation, and, according to Eq. (I2.ip . a ~ 2x + 0(x 2 ) = 2(r/ri). Furthermore, as it 
can be easily appreciated, a = Pm/Pr- Slightly different conventions can be adopted 0. 

Defining as T; the initial integration time (which will be after neutrino decoupling but before matter- 
radiation equality) it will be required that kr- x < 1 for all the modes involved in the calculations. The 
smallness of kr\ measures the excess of the given wavelength with respect to the Hubble radius at the 
time T\ (typically selected well before r cq ). Summarizing, the double expansion employed in setting 
initial conditions of the truncated Einstein-Boltzmann hierarchy can be formally written as 

k k na , „ ,,,2 y/a + 1 . 

a<l, — = - = - 1 ===~kT<l, H=-- . 2.2 

aH H 2 V a + 1 T\ a 

where k = kr\ measures how large the wavelength was, in Hubble units, around equality (note, indeed, 
that r eq = (y/2 - l)n ~ r x /2). To derive the expression of the various adiabatic and non-adiabatic 
modes it is practical to use, as a guiding principle, the evolution of curvature perturbations in the 
limit k/7i ~ kr — > 0. The obtained expression will be used (to leading order in the a-expansion) for 
solving the truncated Einstein-Boltzmann hierarchy to a given order in kr. The two mentioned steps 
can be iterated at wish to obtain the initial conditions with the wanted accuracy either in kr or in a. 

Before plunging into the discussion it is appropriate to recall [7J |TT] that large-scale magnetic fields 
affect directly the Einstein equations since they contribute to the energy density, to the pressure and 
to the anisotropic stress: 

where B 2 = BiB 1 . For notational convenience, and in analogy with what customarily done to treat 
the neutrino anisotropic stress [25J, we will also define 

<),<)' \\j = (p 7 + p 7 ) vV B , ms, t) = 4pb( , x ; t) , (2.4) 

where <7b and ^b are both dimensionless. The magnetic component also enter directly the Boltzmann 
hierarchy. In the truncated version, which is used to set initial conditions the magnetohydrodynamical 
(MHD) approximation, the Lorentz force affects the baryon-lepton-photon evolution. In this respect 
it is useful to bear in mind that 



V-(Jx5) 1 
a 4 p h R h 



(2.5) 



where J ~ V x B/(Att) is the Ohmic current and is the baryon-to-photon ratio which can be 
expressed as 

3 p b 3 n M tJbo 3a Lo h0 1 a 

4 p 7 4it 7 \ 2 R uuo 4it 7 ujmo 1 + a 1 + a 

where, by definition, R 1 = p 7 /pr. The relevance of i?b depends on the specific solution. As it can be 
also appreciated from Eq. (1A.18|) . the contribution of i?b is always suppressed, before equality, since it 



10 In the literature there is sometimes the habit of setting t\ = 1/2 and a cq = 1/4. This choice would imply that 
a(r) = r + r 2 . In the latter case the asymptotic solution during radiation will go as a(r) ~ r. Thus the multiplicative 
factors of the various multipoles will change in comparison with the conventions of the present paper. In the numerical 
code, the evolution of the scale factor is computed numerically as a function of r measured in comoving units of Mpc 
(while k, the wavenumber, will be measured in units of Mpc -1 ). 
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is proportional both to a C 1 and to Wbo/^Mo- However, to get accurate initial conditions it is often 
important to take into account the contribution of to the asymptotic form of the evolution equation 
of the baryon-photon dipole. To leading order in k/Ti the evolution of the curvature perturbations^ 
reads [6j [7] : 

n' = + -^-(4 - \)s bPb + o(£ s \ (2.7) 

Pt+Pt Pt + PtV 3/ \H l ) 

The remaining terms appearing in Eq. (1A.23I) will be, as specified in Eq. (12.71) . at most, 0(k 2 T 2 ) 
except for the peculiar case of the neutrino- velocity mode where they will be, at most, O(kr). Once 
the solution of TZ is known, the obtained result can be used in the following equation: 

Equation (12. 8ft . read from right to left, is just the definition of TZ in terms of the synchronous degrees 
of freedom. However, if TZ is determined from Eq. (12.71) . Eq. (12.81) allows to determine £. The solution 
of (I2.8P can finally be inserted into the Hamiltonian constraint, i.e. Eq. (]A.2jl : 

2k 2 i-Hti = 8nGa 2 [6 s p t + 5 s p B }, (2.9) 

to determine h(k,a). Since both Eqs. (12. 8p and (12.91) are exact, the solution of Eq. (12.91) even allows 
for the determination of h(k, a) to first order in k 2 r 2 . Across the radiation-matter transition, the total 
barotropic index and the total sound speed can be written, respectively, as 

w t = — = ^ — -TV, c st = — = — — — . 2.10 

p t 3(a + 1) apt 3(3a + 4) 

Therefore, Eqs. (12.71) and (12. 8p become, in explicit terms, 

dTZ AV^k) 
~da ~ ~(3a + 4) 2 ' 



; 2 .n) 



d£ 3a + 4 . 3a + 4 ^ , n ^ n . 
da 2a(a + l) s 2a(a + 1) ' v ' ; 



where, for short, the quantity 



£>*(fc) 



ujmo 4 



(2-13) 



has been introduced. In Eq. (12.111) only one CDM-radiation mode has been included as an illustration 
of the procedure and, in this case, 

3 

4Pnad = p c C 2 t<S*( k )i = -5 r (k, t) - 8 c (k, t). (2.14) 

Equations ( 12. 1 If) and ( 12. 12ft can then be integrated once with respect to a; the result is: 

^,a)=7^)-^S (2.15) 
3a + 4 

a) = lZ*(k) - ^^[a(a - 4) + 8( v / aTT - 1)], (2.16) 



11 The evolution equation of the curvature perturbations is swiftly derived, for completeness, in Eq. (|A.23[) of Appendix 
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In the long wavelength limit defined by Eq. (12.21) . the density contrasts following from Eqs. (1A.13j) - 

( ELBD 

S u (k, a) ~ -B^Qsik) + \h(k, a), (2.17) 
5 7 (A;, a) ~ -i? 7 ft B 0) + \h(k, a), (2.18) 



6 b (k,a) ~ -^fieCfc) + ^rr^' ( 2 - 19 ) 
<J c (Ar, a) ~ - -^fieO) + (2.20) 



3 D n /-m , MM) 
4" 

If S*(k) = in Eq. fl2T20l . Eqs. fI^TT|) f|272Q|) will describe the adiabatic mode. If S*(k) ^ in Eq. 
( 12.201) . Eqs. ( 12.17l) -( l2.20l) apply in the case where the non-adiabatic contribution is present either in 
combination with the adiabatic mode or in the absence of the adiabatic component. 

By inserting Eq. (12.161) into Eq. (I2.9P (written in the form given in Eq. (IA.8I) ) a decoupled equation 
for h(k, a) can be obtained with the combined use of Eqs. (I2.17p - (l2.20p : 

eh 3^ + 4 mwg^ 

da 2a(a + l) a + 1 H 2 a v ; 

Equation (12.211) can then be solved with elementary methods and the result is: 

h(k, a) = V^T^a) + K 2 1l*(k) F 2 (a) - n 2 V*(k) JF 3 (ai), (2.22) 

where 



= 2 [ n ( n -4) + 8(Vo + l-l)] 3 
ar 



= [a »- 2^ + to + 16(1-^+1)] 
5a 2 



{32(V^+l-l) + a[a(3a- 26) + 60 v / a + l-16]-60 v / « + lln(a + l)} . 
^3 « = Tr~ 2 • 2-25 

In the absence of magnetized contribution TZ(k, a) is constant, to leading order in a and kr, provided 
S*(k) = 0. This is also true for £(k,a). The situation for h(k,a) is a bit different. Indeed, the 
synchronous coordinate system implies that the gauge freedom is not completely fixed. Thus the 
possible presence of a constant mode in h(k, a) is the reflection of a gauge artifact which can be 
removed by a further gauge transformation which will keep the transformed fluctuation always within 
the synchronous coordinate system [241 1451 147] . In our approach the gauge mode does not show up 
since we integrated directly the evolution equation of curvature perturbations (i.e. Eq. (|2.7|) ) which is 
gauge-invariant and hence free of possible gauge modes. By selecting a given order in a, the asymptotic 
solution can be fed back into all the other equations of the truncated Einstein-Boltzmann hierarchy. 
The system can then be solved to a given order in kr. There is no general way of achieving the second 
step. It will thus be mandatory to solve the system separately for the different modes. 



3 Magnetized adiabatic mode 

The long wavelength solution derived in the previous section implies that, to leading order, 

r) ~ K^k) + 0(k 2 r 2 ), r) ^ Uk, t) ~ ( c (k, r) ~ Cb(*, r). (3.1) 
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Figure 1: The magnetized adiabatic mode for spatially flat models is illustrated. The ACDM param- 
eters are selected to match with the ones of the WMAP three year best fit to the WMAP alone (i.e. 
h = 0.732, uj w = 0.02229, uj c0 = 0.1054 and An = 2.35 x 10~ 9 ). The combination of the WMAP 
data with the other sets of cosmological data (or with the other CMB experiments) leads to best fit 
parameters which are largely compatible with the ones used here only for sake of illustration. 



Thus, to leading order in a and in kr the <Sy introduced in Eqs. fl 1 -2f> and (11.31) are all vanishing. The 
result of the integration of the truncated Einstein-Boltzmann hierarchy of Appendix El can then be 
written, in Fourier space, as: 



h(k,r) 



6 u (k,r) 
6 c (k,r) 



+ 



ARu + 5 K (k) + ^t 4 ^) ~ R Mk)] 



12(AR U + 15) 
1 



8Rl 



IAR V - 75 



36 1 2(2^ + 25X4/^ + 15) 



6(4i^ + 15) 
K*{k) 



2^2 



10(4^ + 15X2^ + 25) 
-R^k) + - 



[RMk)-4a B (k)]^r\ 
k 2 r\ 



^ + a B (k) - ^n B (k) 



-A r n B (fc) + - 



--R,n B (k) + i 
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fcV 
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^a B (k) - ^n B (k) 

Hy{AR v + 27) 
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0, 



w (4^ + 23) 

36 I 4/^ + 15 V ; 



21Z*(k) R^[4a B (k) - R u Q B (k)} 



3(AR U + 15) 



+ 



2/^(4/^ + 15) 



fcV. 



(3.7) 
(3.8) 

(3.9) 
(3.10) 
(3.11) 



The adiabatic mode is often assigned in terms of a putative C{k) which is related to H*{k) as lZ*{k) = 
—2C(k). To leading order in kr the initial conditions are adiabatic, however, higher orders in kr 
introduce a perturbative mismatch between the relevant density contrasts. Higher orders in a (at 
a given order in kr) can be easily included following the expressions of the previous section. Since 
the Hamiltonian constraint has to be satisfied at the onset of the integration, the contribution of 
the magnetic field in the initial data acquires a definite sign. The relative correlation between the 
magnetic field and the adiabatic mode is fixed. In Fig. [TJ the salient features of the adiabatic mode 
are illustrated and have been more thoroughly discussed in [TTJ. For sake of comparison with the 
other sets of initial conditions, the adiabatic power spectra are assigned according to the standard 
convention, i.e. Vn(k) = An(k / 'k p ) ns ~ l , where k p = 0.002 Mpc -1 is the pivot scale, i.e. the scale at 
which V n (k p ) = A n - The three year best fit to the WMAP data alone implies [H[2] A n = 2.35 x 10~ 9 . 
According to the legends of each plot, in Fig. [T]the full line denotes the corresponding angular power 
spectra in the absence of magnetic contribution and for WMAP three year best fit. 

The magnetic fields are included in the simplest realization of what has been called in [TU] mACDM 
model, i.e. the magnetized ACDM model. In this case the power spectra of Q B and a B are assigned, 
respectively, as Vn{k) = VL\^(n B )(k/k h ) 2 ^-^ and as V a {k) = U 2 BL g (n B ) (k / ^A; L ) 2(nB_1) , where U BL = 
Bl/ '(87ra 4 p 7 ). The functions ^(^b) and Q{n B ) can be determined from the appropriate convolutions 
[TT) once the two-point function of the magnetic fields is known in Fourier space, i.e. 



2tt 



(Bi{k)Bj(p)) = —P l3 (k)V B (k)5^\k + p), P^k) 



(3.12) 



The quantity is the amplitude of the magnetic field regularized over a comoving scale L ~ k^ 1 . In 
this context f2 B L measures (up to Euler Gamma functions [TTJ) the amplitude of the magnetic power 
spectrum V B (k) at the magnetic pivot scale k-^. From Fig. [TJ the TT spectra exhibit a distortion of 
the second peak which is correlated with the increase of the first and third acoustic peaks. As the 
values of the magnetic spectral index (i.e. n B ) increases the distortion becomes more pronounced. 
As the values of the magnetic field augments beyond 20 nG the second peak practically disappears 
and it is replaced by a sort of hump. This result is also expected on analytical ground [6j O [8] . In 
Fig. [TJ excessive values of the magnetic fields have been adopted just to emphasize more visually the 
modifications of the shape which are typically related to this mode. This aspect can be understood 
by comparing, in Fig. [ljthe dashed (and dot-dashed) lines with the full lines (representing the three 
year best fit to the WMAP alone). In the bottom-right plot the different curves denote three different 
values of the magnetic fields for fixed magnetic spectral index (i.e. n B = 1.2) and fixed adiabatic 
spectral index (i.e. n s = 0.958). 
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The magnetic pivot scale is taken (in Fig. [T] and in all the other plots of this paper) to be 
1 Mpc -1 . The EE and TE correlations are also modified by the presence of large-scale magnetic 
fields. The physical rationale for this occurrence can be understood already to lowest order in the 
tight-coupling approximation: since the dipole is affected both by a B and fl B the TE and EE angular 
power spectra are modified since the polarization (to first-order in the tight-coupling expansion) is 
proportional to the zeroth-order dipole which depends, in the MHD treatment, on the Lorentz force 
and, ultimately, on the Ohmic currrent. This aspect has been already pointed out, within a much less 
accurate semi-analytical approach, in [6]. 



4 Magnetized C DM- radiation mode 

In the case of the magnetized CDM-radiation mode the entropy fluctuation arises in the CDM-radiation 
system. Following the notations of Eq. (11. 5p we will denote by S*(k) the common amplitude, in Fourier 
space, of the two relevant entropy fluctuations, i.e. 



S*(k) = Scj(k) = S cv (k). 



(4.1) 



In principle we might have the physical situation where S cl {k) ^ S cu (k). In the latter case, however, 
Eq. (11. 5ft implies that S vl {k) ^ 0. This means we would not be dealing with a pure CDM-radiation 
mode but rather with a CDM-radiation mode supplemented by a neutrino-entropy mode (see also Eq. 
(11.81) ). For the magnetized adiabatic mode the a-expansion defined in Eq. (I2.2p only provides useful 
corrections to the leading order (constant) solution for £(/c, r). For the CDM-radiation solution £ is 
not constant but it rather vanishes as a. Consequently, for a C 1 and fcr<l, the asymptotic form 
of the solution is: 



h(k, t) = h(k, t) = 
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where, as established in Eq. (12. II) . x = t/ti and T>*(k) has been already defined in Eq. (12.131) . The 
function h(k, r) can be expanded for a <C 1: 



h{k,r) 



Va + T 



Q 



16 + 



2(a 2 -4a -8) 



Va + 1 



£>*(fc) 



a 



a 2 m,(fc) + C(« 3 ). 



(4.13) 



The CDM-radiation power spectra are assigned as Vs(k) = As{k j l k^) n ^ x with the same pivot scale, 
i.e. k p = 0.002 Mpc^ 1 adopted in the adiabatic case. In Fig. [2] the spectral index of each plot (i.e. 
n c ) is specified in the corresponding legends^]. To ease the comparison between Figs. [T] and [2] the 
non-adiabatic amplitude is fixed to the same value of the adiabatic one, i.e. As = -Mi — 2.35 x 10~ 9 . 
For n c > 1 a sort of hump arises in the Doppler region. In Fig. [2] (top-left plot) the TT correlations 
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L B c 

B L = 10 nG, n B = 1.5, n -2 
B, = 10 nG, n = 1.5, n =3 
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Figure 2: The magnetized CDM-radiation mode is illustrated for spatially flat models. The cosmolog- 
ical parameters, for sake of comparison, are fixed to the same values of Fig. [TJ 

are reported in the case n c = 2. The full line denotes in this plot the case = while the dashed and 
the dot dashed lines correspond to the case = 10 nG with two different magnetic spectral indices 
(as indicated in the legends). The top- right and bottom-left plots of Fig. [2] the EE and TE angular 
power spectra are illustrated with the same notations. In the bottom-right plot the TT correlations 

12 In the title of the plots the value of the optical depth has been also specified (i.e. r = 0.089 as implied by the 
WMAP three year best fit). This notation is standard and we did not deviate from it. However, it should be remarked 
that t is consistently used in our script to denote the conformal time coordinate. This remark avoids any unwanted 
confusion. 
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are reported for three different values of the non-adiabatic spectral index (i.e. n c = 1, 2, 3) when 
the parameters of the magnetized bckground are fixed. From the latter plot it is clear that, in the 
scale- invariant limit the hump is not present. 

After careful numerical scrutiny of the magnetized CDM-radiation mode the following general 
aspects can be highlighted: 

• the inclusion of the magnetic field diminishes the height of the hump of the TT angular power 
spectra and also suppresses the TE and EE correlations; 

• slight changes in the magnetic spectral index (e.g. from 1.2 to 2 in Fig. [2]) affect quantitatively 
all the CMB observables; 

• the phases of the EE angular power spectra (and, to a lesser degree, of the TE angular power 
spectra) are slightly shifted as £ > 500. 

The TE correlation starts positive at large scales (see Fig. [2]bottom-left plot). This is in sharp contrast 
with what is experimentally observed (i.e. the TE correlations exhibit an anticorrelation peak [TJEJEj). 
The presence of the magnetic field can indeed reduce the TE amplitude but cannot change the sign of 
the first peak by making it negative. 

5 Magnetized baryon-radiation mode 

If "C ^co ( as h happens in the ACDM paradigm) the baryon-radiation mode is, grossly speaking, 
less important than the CDM-radiation mode. The baryon-radiation mode touches upon the baryon- 
lepton fluid (which is also affected by the Ohmic current). Interesting effects on the phases of the TE 
and EE angular power spectra are then potentially expected. From Eq. (11.61) we will denote by 5*(/c) 
the common amplitude, in Fourier space, of the two relevant entropy fluctuations, i.e. 



For the reasons explained in the previous Section, S^(k) ^ S^ v (k) also demands S vl {k) ^ and this 
would imply that the baryon-radiation mode is supplemented by a neutrino-entropy mode. To make a 
distinction with respect to the CDM-radiation mode, the Fourier amplitude of the entropy fluctuation 
has been denoted, in the baryon-radiation case, as S*(k). Direct integration of the truncated Einstein- 
Boltzmann hierarchy summarized in Appendix [A] leads to the following initial conditions: 



S*(k) = S h ^(k) = S hu {k). 



(5.1) 




(5.2) 



(5.3) 



5 7 (k,r) = -h(k,r) -_R 7 fi B (&), 
6 u (k,r) = h(k,T)-R,Q B (k), 



(5.4) 



(5.5) 



6 c (k,r) = ^fl-^R,Q B (k), 



(5.6) 
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Figure 3: The magnetized baryon-radiation mode is illustrated for spatially flat models. To ease the 
comparison, the cosmological parameters are fixed to the same values of Figs. Q] and EJ 



6 h (k,r) 
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R 
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where, following the analog notation used in Eq. (12.131) for the CDM-radiation mode, 

^bO ^bO 



Wmo 4 



(5.7) 

(5.8) 

(5.9) 
(5.10) 
(5.11) 

(5.12) 
(5.13) 



^MO ^cO + ^bO 

The power spectrum of the baryon-radiation mode is assigned as Vg{k) = Ag(k / k x> ) flh ~ 1 where 
is the corresponding spectral index. Comparing Eqs. fl5.2l) - fl5.12l) with Eqs. fl4.2p - fl4.12p . there are 
formal analogies but also some differences. In the case of Eqs. fl5.2p - fl5.12p the baryons are directly 
affected and, therefore, 
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• as anticipated, for the fiducial set of parameters of the ACDM model the amplitude of the 
TT, EE and TE angular power spectra will be comparatively smaller than in the case of the 
CDM-radiaation mode; 

• a stronger effect on the phases is expected when the magnetic field is consistently included in 
the calculation. 

In Fig. [3] the TT, EE and TE angular power spectra are illustrated when the amplitude is selected in 
such a way that A$ = An- The comparison of Fig. [3] and [2] suggests that, indeed the magnetized CDM- 
radiation mode is very similar to the magnetized baryon-radiation mode. Also in the case of the baryon 
radiation-mode the hump arises for > 1 The solution of the Hamiltonian and of the momentum 
constraints implies a definite relative sign of the magnetic and of the entropic contributions. As a 
consequence the amplitude of the angular power spectra diminishes with respect to the case = 0. 
The phases of oscillations of the TE and EE angular power spectra (see top-right and bottom-left plots 
in Fig. [3]) are shifted as the magnetic field strength increases. 



6 Magnetized neutrino-entropy mode 

In the case of the magnetized neutrino-entropy mode the non-adiabatic fluctuation belongs to the 
neutrino sector, i.e. 

S m (k) = l[6 1 (k,T)-6 v (kr)]. (6.1) 

The lowest order form of the solution can be obtained by setting £(k,r) = h{k,r) = 0. This aspect 
already highlights the difference of the neutrino-entropy mode in comparison with both, the previous 
non-adiabatic modes and with the adiabatic solution. In the case of the adiabatic solution, to lowest 
order in the expansion of Eq. (12.2ft . the curvature perturbations are constant in time, i.e. 



K(k, a) ~ K*(k) + 0(k 2 r 2 ) + 0(a). 



(6.2) 



For the CDM and baryon-radiation modes the curvature perturbations vanish as a, to lowest order in 
fcr, i.e. K oc a + 0(k 2 T 2 ). In the case of the neutrino-entropy mode the behaviour is intermediate since 
TZ(k, a) ~ 0(k 2 r 2 ). The specific form of the neutrino-entropy mode is then obtained by integration of 
the truncated Einstein-Boltzmann hierarchy reported in Appendix [A] The final result is: 
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(6.7) 
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Figure 4: The angular power spectra obtained in the case of the neutrino-entropy mode are illustrated 
for the same set of cosmological parameters used in Figs. [TJ El and El 
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(6.11) 
(6.12) 

\(n„-l) 



In the case of the neutrino-entropy mode the power spectra are assigned as Vg = Ag(k/k f 
where n v is the corresponding spectral index. In Fig. [5] the numerical integration of the angular power 
spectra is reported for Ag = An = 2.35 x 10~ 9 . Figure H] also shows two interesting features: 

• the hump disappears not only in the scale-invariant limit (i.e. n u = 1) but also for n u > 1; in 
the latter cases, however, the amplitude increases; 

• as in the other two cases the Hamiltonian and momentum constraints imply that the inclusion 
of the magnetic field always tend to diminish the heights of the acoustic oscillations. 



As it is typical also of the other modes the TE angular power spectrum always exhibits a correlation 
(rather than anticorrelation) at large-scales. The phase shifts in the polarization observables do not 
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arise for the neutrino-entropy mode: the magnetic and the entropic contributions do not interfere 
directly since the neutrinos, unlike the baryons and electrons, only feel the magnetic field through the 
constraints and through the ansiotropic stress. 

For sake of completeness, the solution of the truncated Einstein-Boltzmann hierarchy will be derived 
for the neutrino-velocity mode. The neutrino-velocity mode is rather peculiar since it is, strictly 
speaking, not an entropic mode. It is derived by requiring that the initial fluctuation resides in 
the velocity of the neutrinos. The momentum constraint stemming from the (Oi) component of the 
perturbed Einstein equations (i.e. Eq. (IA.9I) ) then demands that also the photon-baryon velocity is 
perturbed at early times so that, overall, the momentum constraint is satisfied at the onset of the 
numerical integration. From this observation the leading order solution implies that h ~ and £ ~ 
and, therefore, 

i2„0„ + -Ry(l+flb)ftyb->O, (6.13) 

at early times. Of course, both 9 U and # 7 b are both non-vanishing. The neutrino-velocity mode can 
be initialized by setting: 

h(k,T)=-ls*(k)^^kT lX 2 , (6.14) 

iS?^ (6 - i5) 

*r(*. T ) = ~S m (k)hr - SM l) ^kr lX 2 - R,n B , (6.16) 

6 rL/ xt 7 Co>mo 

S v (k, t) = -^-S*(k)kT - S^(k)^-^-kTix 2 - R 7 Q B (6-17) 

3 Ky o>mo 

5 c (k ) r) = --Mk)^^kr 1 x 2 (6.18) 

4 Ky Wmo 

5 h (k,r) = ^(k)kr - -Mk) K{R 2 2 + ^ fcV, (6.19) 
riy 4 _« 7 u; M0 

MM = -^(k) + 3 J^-)^(k)kx - !§(^Y l Uk)kx 2 
R~f 2 Vc^mo/ Ry 4 Vwmo/ 
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k l r (6.21) 
9 C = 0, (6.22) 
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Ry 45* (fc) 



^, r) = -_2 ffB + _^ T , (6 . 23 ) 

where, as established in Eq. (|2.ip . x = t/t\. The quantity S*(k) arising in the neutrino-velocity 
mode is physically different from the analog variable appearing in the solution of the neutrino-entropy 
mode. Here it measures the spectrum of 9 U (or 7 b). While the usual non-adiabatic modes can be 
described also in different gauges, the neutrino-velocity mode might even diverge if we move from the 
synchronous gauge. In the longitudinal gauge the scalar fluctuations of the geometry are defined, in 
the spatially flat case, as 5 s g^ = 2a 2 (p and as S s gij = 2a 2 ip5ij. The relation between the longitudinal 
and the synchronous degrees of freedom is given by 

[(h + QO" + n(h + QO'} C| M(ft + 6gy (g24) 



2k 2 ' 2k 2 
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From Eqs. (16.31) and (16. 4p the values of <p and ip can then be obtained to leading order in kr and they 
are: 



<j>{k) 



3(4i^ + 15) 



S*(k) - ^-a B (k) + -Mk) 



(6.25) 



The same exercise can be performed in the case of the neutrino- velocity mode. Thus, inserting Eqs. 
(IBTTIj) and flBTTS]) into Eq. (KT23jl we do get that 



<f>(k,r) 



AR u S,{k) 



if)(k,r) = -(f>{k,r). 



(6.26) 



(4^ + 5)1^1' 

Consequently, the neutrino-velocity mode is singular in the longitudinal gauge. In the same gauge 
the neutrino-entropy mode is not singular. It is therefore mandatory to use the synchronous frame 
to describe the neutrino-velocity mode. The unphysical nature of this divergence is clear since the 
gauge-invaraint curvature perturbations, i.e. 7Z, is regular for both modes. 



7 Mixtures of magnetized initial conditions 

Assuming a precise knowledge of the thermal history of the Universe from the Planck energy scale 
down to photon decoupling, then, definite predictions on the CMB observables can be made. To 
make the predictions more definite, it is wise to assume that the field content during inflation is 
rather limited. If we just demand a single adiabatic mode, as in the case of the ACDM paradigm, 
then a single inflaton field must be present. This simple model is constrained by experimental data. 
More complicated models will be more constrained as soon as more accurate data will be available [9]. 
When magnetic fields are consistently included in the pre-equality physics, the issue of initial conditions 
becomes particularly important and it has been recently shown that different thermal histories can 
imply slightly different amplitudes of the curvature perturbations induced by a magnetized background 

gg. 

In a pragmatic perspective it is plausible to start confronting with observations the simplest model. 
In the case of [10], the mACDM scenario is the simplest option since a magnetized adiabatic mode is 
assumed and the magnetic fields introduce only two new parameters in comparison with the ACDM 
paradigm. It is however equally plausible to bear in mind also different possibilities. 

The most general initial conditions of the magnetized CMB anisotropies will be given, presumably, 
as a mixture of various modes. As previously observed the TT angular power spectra will be typically 
distorted and shifted upwards if only the magnetized adiabatic mode is assumed in the pre-equality 
plasma. The magnetized CDM-radiation mode tends to reduce the amplitude of the TT power spectra 
in comparison with the case when the magnetic fields are absent. Since the magnetized adiabatic mode 
and the CDM-radiation mode are both solutions of the truncated Einstein-Boltzmann hierarchy, also 
their combination will lead to a viable initial condition. The same reasoning can be applied also in 
more complicated cases, like the one where, on top of the magnetized adiabatic modes there are two or 
three non-adiabatic modes. Generally speaking, the relative amplitudes of the different contributions 
must be gauged in such a way that: 

• the acoustic peaks are correctly reproduced in the TT angular power spectrum; 

• the anticorrelation peak in the TE spectra is not erased by the addition of the non-adiabatic 
component; 
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Mixed initial conditions with n B =1 .5, B L =1 nG 
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Figure 5: The TT angular power spectra in three different cases: AD (magnetized adiabatic mode), 
AD + NAC +NAB (mixture of the magnetized adiabatic mode supplemented by the magnetized 
CDM-radiation and baryon-radiation modes), AD + NAC + NAB +NE (same as in the previous case 
but with the addition of a neutrino-entropy component). The spectral parameters are fixed as in Eqs. 
(TTTj). CL2D and (Q. 



• the EE angular power spectra should not be too different, at large multipoles, from the extrap- 
olated three year best fit to the WMAP data alone. 

The third requirement can be phrased in terms of any best fit (also including other data sets, e.g. 
[T21 [T5] and [TU US])- This will not change our conclusions since, as far as CMB observables are 
concerned, different best fits lead to central values that are not crucially different for the present 
(illustrative) purposes. In Fig. [5] the TT correlations are reported for different sets of initial conditions. 
In the top-left plot the full range of multipoles is discussed and in the remaining three plots the low, 
intermediate and large i regions are specifically scrutinized. As the legends specify the full thick line 
denoted the three year best fit to the WMAP data alone. The cosmological parameters (i.e. c^bO; ^co> 
n s , r and An) are hence fixed to the values already assumed for the previous plots. The light curve 
denotes the magnetized adiabatic mode in the case = 10 nG and n-s = 1.5. Inspection of the plots 
reported in Fig. [5] suggests that the case labeled by AD (i.e. a single magnetized adiabatic mode) is 
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incompatible with the experimental data. This (expected) conclusion can be partially evaded if the 
magnetized adiabatic mode is combined with a mixture of initial conditions. Two different cases will 
then be considered: 

• AD + NAC+NAB denotes the case where, on top of the magnetized adiabatic mode (i.e. AD), 
the initial conditions include also two other non-adiabatic contributions which are chosen to be 
the magnetized CDM-radiation (i.e. NAC) and the magnetized baryon-radiation (i.e. NAB) 
modes derived, respectively, in Eqs. (14.2p - (14.12j) and in Eqs. f l5.2p - fl5.12p ; 

• AD+ NAC+ NAB+ NE denotes the case where, on top of the three contributions mentioned in 
the previous point there is also a neutrino-entropy mode (i.e. NE) which has been derived in 
Eqs. fl6\3|) - fl67[2|) . 

In the case AD + NAC + NAB the spectral parameters are fixed, according to the present notations 
as: 

As = 0.1 A n , n c = 1.5, Ag = 0.3 A n , n h = 1.5. (7.1) 

For the case AD+ NAC+NAB+NE, the spectral parameters of the neutrino-entropy mode are fixed 

as 

Ag = 0.025^, n v = 1.4, (7.2) 

while the spectral parameters of the other two non-adiabatic modes are chosen as in Eq. (17. II) . For 
both examples summarized in Eqs. (17. ip and (17. 2p the remaining spectral parameters are selected as: 

J B L = 10nG, n B = 1.5, A n = 2.35 x 1(T 9 , (7.3) 

The values of the other cosmological parameters is the one suggested by the three year best fit to the 
WMAP data alone. According to Fig. [51 the excessive increase (and distortion) of the acoustic peaks 
induced by the bare magnetized adiabatic mode (i.e. AD) is compensated when the initial conditions 
are of the type AD + NAC + NAB. The NAC contribution, as discussed in Fig. [2] tends to compensate 
the increase of the peaks while the NAB contribution acts on the phases (see also Fig. [3]). As it is clear 
from the bottom- right plot of Fig. [5] the AD+NAC+NAB combination is even too effective in lowering 
the peaks. The inclusion of the neutrino-entropy mode moderates this trend by adding extra-power in 
the acoustic region. This effect arises since the neutrino-entropy mode has a different phase structure 
in comparison with the non-adiabatic mode of the CDM and baryon sector. From the different ranges 
of multipoles where the TT power spectra have been computed two opposite remarks seem to emerge: 

• on the one hand the mixing of different magnetized modes in the pre-equality initial conditions 
seems to shade the effect of a (pretty strong) magnetic field; 

• on the other hand the reported WMAP error bars also suggest that the different cases can still 
be distinguished. 

In the light of these remarks it is useful to consider the plots reported in Fig. [6] where the TE and 
EE angular power spectra are illustrated. The sub-dominant non adiabatic components imply that 
the leading effect in the TE anticorrelation peak is determined by the magnetized adiabatic mode. 
At the same time the sharp increase of the peaks in the TE and EE spectra (for high multipoles) is 
partially compensated by the non-adiabatic components. For £ > 1000 the distortions induced by the 
magnetized adiabatic mode can be reduced but not eliminated. 
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Mixed initial conditions with n B =1.5, B L =10 nG 
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Figure 6: The TE and EE angular power spectra are illustrated for the same mixtures of initial 
conditions reported in Fig. 



The reported examples are, in our opinion, rather effective for motivating the approach pursued 
in the present research program which is directly linked to the possibilities of achieving accurate 
determinations of the temperature and polarization observables with the forthcoming Planck explorer 
mission [!9j. A systematic approach in this direction is a necessary step for any quantitative strategy 
of parameter extraction involving large-scale magnetic fields. 



8 Concluding discussion 

More than fifty years ago large-scale magnetic fields entered the scene because of the interplay with 
cosmic ray physics [49]. Today cosmic ray physics is again confronted with the relevance of large-scale 
magnetic fields to infer more fundamental conclusions on the nature of the observed fluxes of ultra- 
high energy cosmic rays [50]. In the intervening lapse of time radio-astronomy has established that 
the largest gravitationally bound systems are all magnetized to various degrees. Still the question 
remains: where do the observed fields come from? 

The standard paradigm used to interpret the cosmological data, i.e. the ACDM model, does not 
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contemplate the presence of a magnetized contribution very much in the same way as, in the pivotal 
ACDM mode tensors are absent. The increase in the quality of astrophysical observations will allow, 
in the near future, to include further parameters to the basic ACDM scenario. In [10] we formulated 
a plausible extension of the ACDM scenario, i.e the magnetized ACDM model (mACDM). In the 
minimal mACDM model the CMB initial conditions are simply dictated by the magnetized adiabatic 
mode where no entropic modes are allowed. In the present paper the minimal mACDM model has 
been extended to the case of magnetized initial conditions with entropic components. The spirit of 
the present investigation is, at once, modest and conservative: we just want to bring the treatment of 
magnetized CMB anisotropies to the same standards which are typical of the cases where large-scale 
magnetic fields are absent. The problem with this plan is that specific numerical tool have to be 
developed. Unlike more standard additions to the ACDM paradigm (which are customarily included 
in the standard codes) large-scale magnetic fields demand dedicated numerical approaches which, so 
far, did not exist. The proof of this statement is that, before the present analysis, we did not know, for 
instance, the effects of large-scale magnetic fields on the various non-adiabatic modes contemplated in 
the standard CMB initial conditions. 

In the present paper the initial conditions of magnetized CMB anisotropies have been derived in 
a generalized perspective. The analysis has been conducted both analytically on the truncated form 
of the Einstein-Boltzmann hierarchy. The CMB observables have then been numerically computed 
for each single set of initial conditions to emphasize the relevant shape effects. The addition of the 
magnetic field is not able to change the sign of the first peak of the TE correlations so that we can 
safely conclude that a single magnetized non-adiabatic mode is not compatible with the experimental 
observations. The same sharp statement cannot be made, however, if the initial conditions of the 
truncated Einstein-Boltzmann hierarchy are given in terms of a superposition of modes where the 
dominant amplitude is provided by the magnetized adiabatic mode. The examples provided in this 
paper suggest that a thorough analysis of the TT, TE and EE angular power spectra allow, in principle, 
to distinguish the signatures of different magnetized initial conditions. In the near future we plan to 
pursue our program by sharpening those theoretical tools that might allow for a direct observational 
test of the effects of large-scale magnetic fields on CMB anisotropies. While we are aware of the 
theoretical challenges associated with this effort, the potential relevance of the whole strategy justifies 
our endeavors. 
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A Truncated Einstein-Boltzmann hierarchies 



The evolution equations of the truncated Einstein-Boltzmann hierarchy will be illustrated in the tightly 
coupled regime which is the relevant one to set initial conditions. For reasons of convenience the 
language of the synchronous gauge will be adopted. Similar equations discussed either in different 
gauges or in dedicated gauge-invariant formalisms can be usefully found elsewhere [HJ [T], [8] . The 
synchronous coordinate system is defined from the perturbed form of the metric tensor in terms of 
two variables denoted, in Fourier space, as h(k,r) and r): 



$s9ij (hi 



a 



kikjh(k, t) + 6( kjtj — — )£(k,T) 



(a.i; 



The bulk velocity of the plasma 9^ coincides with the velocity of the baryons which are coupled 
to the electrons through Coulomb scattering. Also the photons are coupled to the baryons because 
of Thompson scattering. Deep in the radiation-dominated epoch, when the initial conditions of the 
Boltzmann hierarchy are ideally set, the baryon velocity field 9^ and the photon velocity field (i.e. # 7 ) 
are very well synchronized because the Thompson rate is much larger, at the corresponding epoch, 
than the Hubble rate. To zeroth-order in the tight-coupling expansion, therefore, 9 1 ~ #b = 9^h where 
# 7 b denotes the common value of the baryon-photon velocity. At later times, well after equality, the 
primeval plasma recombine and the baryon velocity will no longer equal the photon velocity. This 
effect is taken into account in the numerical integration, however, at early time, the zeroth-order in 
the tight-coupling limit is an excellent approximation and will be the one used to set, consistently the 
initial conditions of the Einstein-Boltzmann hierarchy. Higher-orders in the tight-coupling expansion 
can be used also for semi-analytical estimates [HIE]- 

In the synchronous gauge the perturbed Einstein equations for the system of magnetized pertur- 
bations read [TT] : 

2k 2 £ - Hti = 87iGa 2 (5 sPt + 5 s p B ), (A.2) 
k 2 £' = -AnGa 2 (p t + p t )9 t , (A.3) 
h" + 2Uh! - 2k 2 £ = 24nGa 2 (5 s p t + 5 s p B ), (A.4) 
(h + 60" + 2H(h + 60' - 2k 2 £ = 24nGa 2 [(p u + p v )o v + (p 7 + Pl )a B ], (A.5) 

where the following global variables have been defined: 

(Pt + Pt)9t = XX^a + Pa)#a, 6 a fh = ^Pa, *aPt = S sPa. = W a 5 a p a . (A. 6) 



The sums appearing in Eq. flA.6j) extends over the four species of the plasma (i.e. photons, neutrinos, 



baryons and CDM particles) and w a is the barotropic index of each species. Equations flA.21) and 
flA.3j) are, respectively, the Hamiltonian and momentum constraints coming from the (00) and (Oi) 
components of the perturbed Einstein equations. By using the Friedmann-Lemaitre equations in the 
form: 

m 2 = 8nGa 2 p t} H 2 -H' = 4vrGa 2 (p t + p t ) (A.7) 



13 Within the notations of the present paper 8x = V ■ vx is the divergence of the peculiar velocity of a given species 
labeled by X. In the text we will simply talk about "velocity field" but it is understood that this terminology denotes 
the three-divergence of the velocity field. 



23 



Eqs. (1A.2I) and (1A.3j) can also be written in more explicit terms as 



2k 2 i - Hh! = m 2 \ n R (A y 5 7 + r„5 v ) + R y n R n B + n 



M 



<^c0 



-2H 2 



\Cg>M0/ 



V+(^>J}, (A.8) 

(A.9) 



With analog algebra we will have that Eqs. ( 1A.4I) and ( 1A.5I) can be recast in the following form: 



h" + 27tfi' - 2k 2 £ = 3ft 2 fi R (i? 7 5 7 + #A + i? 7 B 



7"7 

(/i + 60" + 2H(h + 60' - 2k 2 £ = l2H 2 tt K [R u o u + R 7 a B 



where 



1 



a 



(A.10) 
(A.ll) 

(A.12) 

Note that, in the momentum constraint (i.e. Eq. (1A.9|) ) we assumed already the tight-coupling regime 
insofar as # 7 ~ 9^ = 6^. To solve consistently for the initial conditions of the Einstein-Boltzmann 
hierarchy we have also to specify the evolution equations for the monopoles and the dipoles of the 
various species. The evolution equations of the monopoles can be written, respectively, as: 



M 



4 2 
# = -3^+3*', 



Si 



Si 



SI 



4 2 
-367b + 3^ 

b! 



-9 C + 



h! 



The evolution equations for the dipoles are given, instead, as 

k 2 



0' 



-S v - k a u , 



HR h 

6 * + ^RbTT 7b 
e c + w c = 0. 



-&v + 



4(1 + 7 4(1 + i? b ) 



(A.13) 
(A.14) 
(A.15) 
(A.16) 

(A.17) 

(A.18) 
(A.19) 



In Eq. ftKtfb there appears also a v and, consequently, the corresponding evolution equation is given 
by: 

(A.20) 



4 a 3 i-t 



2 4 
—ti - 
15 5^ ' 



15 10 

where is the octupole of the neutrino phase-space distribution. 

Recalling that the fluctuation of the total pressure can be expressed as the sum of the adiabatic 
contribution and the non-adiabatic contribution (see Eq. fll.ip ). the generalized evolution equation 
for the curvature perturbations can be obtained. This equation is used, in Section [2] to get the first 
iteration of the large-scale solution which will allow to determine the various modes. By taking the 
difference between Eqs. (1A.4|) and (1A.5j) we do get that 



+ 2Hg = 4nGa 2 [(p u + p v )a v + (p 7 + p 7 )cx B - S s p t - 5 s p B }. 



(A.21) 
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By now summing up Eq. (1A.21I) and Eq. (1A.8|) (appropriately multiplied by c^) the following equation 
can be obtained: 



e + H(2 + 3c s 2 t )£' = ^c 2 st (h' + QO - k 2 cli 



+4nGa 2 



(Pv + Pv)o v + (p 7 + P 7 )o"b + c st - - )5 s p B - S s p nad 



But now recall that the curvature perturbations in the synchronous gauge simply become: 



n 2 -w 



n , = l_p±_ 

3 pt + Pt 



^ + (2 + 3c s 2 t )£' 



(A.22) 



(A.23) 



where the second relation can be obtained from the first time derivative of 1Z since, by definition, 
c st = Pt/ Pt- Finally, combining Eqs. flA.221) and flA.231) . we do get the wanted equation, i.e. 



n' 



n . n 

OsPnad + 



Pt + Pt 



Pt + Pt 

,2 
st 



H 2 {h' + 6Q 

MGa\ Pt + Pt ) AirGa 2 (p t + p t ) 



+ 



S s Pb 

7-L[{Pv + Pv)Ov + (P 7 + P-yW] 

(Pt + Pt) 



(A.24) 



The terms appearing in the second line of Eq. ( 1 A. 2411 are of higher order in k/7i while the terms 
appearing in the first line are responsible for the leading order solution in the limit fcr C 1. Equation 
flA.24j) will be used in Section [2] to obtain the the evolution of the various modes in the limit of 
wavelengths larger than the Hubble radius. 
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